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We have investigated spin singlet Mott states of spin-one bosons with antiferromagnetic interac- 
tions. These spin singlet states do not break rotational symmetry and exhibit remarkably different 
macroscopic properties compared with nematic Mott states of spin-one bosons. We demonstrate 
that the dynamics of spin singlet Mott states is fully characterized by even- or odd-class quantum 
dimer models. The difference between spin singlet Mott states for even and odd numbers of atoms 
per site can be attributed to a selection rule in the low energy sectors of on-site Hilbert spaces; 
alternatively, it can also be attributed to an effect of Berry's phases on bosonic Mott states. 

We also discuss evidence for spin singlet quantum condensate of spin-one atoms. Our main find- 
ing is that in a projected spin singlet Hilbert space, the low energy physics of spin-one bosons is 
equivalent to that of a Bose-Hubbard model for spinless bosons interacting via Ising gauge fields. 
The other major finding is spin-charge separation in some one- dimensional Mott states. We propose 
charge-e spin singlet superfluid for an odd number of atoms per lattice site and charge- 2e spin singlet 
superfiuid for an even number of atoms per lattice site in one- dimensional lattices. All discussions 
in this article are limited to integer numbers of bosons per site. 
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I. INTRODUCTION 



Spin correlated condensates of spin-one atoms have at- 
tracted considerable interest since 1998 [1,2]. Two-body 
scattering between atoms leads to either antiferromag- 
netic (polar) or ferromagnetic condensates [3-5]. For 
spin-one bosons with antiferromagnetic interactions, the 
condensates break both the rotational symmetry and the 
[/(I) symmetry associated with the phases of the con- 
densates. There are two branches of spin-wave excita- 
tions and one gapless phason mode. For a condensate of 
a finite size, the exact ground state however is a spin sin- 
glet (assuming the total number of atoms is even ) which 
on the other hand is extremely sensitive to external per- 
turbations [5,6]. In rotating traps, two-body scatterings 
furthermore result in spin correlated fractional quantum 
Hall states [7,8]. 

In lattices where the hopping is small, bosonic par- 
ticles can be localized because of repulsive interactions 
and ground states in this limit are Mott insulators [9]. 
Mott-insulating states of cold atoms have been recently 
observed in optical lattices [10]. Spin correlated Mott- 
insulating states in high dimensional optical lattices have 
also been investigated theoretically and reported in [11]. 
Both nematic Mott insulators which break the rotational 
symmetry and spin singlet Mott insulating states were 



found for high dimensional lattices in certain parameter 
regime. Other spin correlated states which break both 
translational and rotational symmetries were proposed 
in [12]. Effects of spin correlations on Mott insulator- 
supcrfluid transitions were suggested and remain to be 
fully understood [13]. Very recently, detailed analysises 
of microscopic wave functions and phase transitions be- 
tween spin singlet Mott insulators and nematic ones have 
been carried out [14,15]. In the context of antiferromag- 
nets, the issue of spin nematic states was raised and ad- 
dressed previously [16-18]. 

To distinguish different correlated states of spin-one 
bosons in a lattice, it is useful to introduce the following 
order parameters in terms of creation and annihilation 
operators V'l, V'a respectively (a = x,y, z); 
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In Eq.l, we have employed the following creation opera- 
tors 
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where V'to^i '"^f = ±1,0, are creation operators of spin- 
one bosons with Sz = mp- 

In terms of these order parameters, polar condensates 

(pBEC), ncmatic Mott insulators (NMI) and spin singlet 
Mott insulating states (SSMI) have the following proper- 
ties 

pBEC :Oi^0,Olp^0; 

NMI :Oi = 0,Olf,^0; 

SSMI:e»i =0,e»2^ = 0. (3) 

Finally, for all states of cold atoms with antiferromag- 
netic interactions which interest us in this paper, the ex- 
pectation value of spin operators is zero 

(S) = 0, S« = -ie«^^V^V7; a, 7 = x, y, x. (4) 

Spin correlations in Mott-insulating states also depend 

on the cvcn-odd parity of numbers of atoms per site. 
While spin correlated Mott states for even numbers of 
atoms per site have been studied and understood, states 
for odd numbers of atoms per site remain to be fully un- 
derstood. In one-dimensional optical lattices, it was ar- 
gued that spin singlet Mott states for an odd number of 
atoms per site should be dimerized-valence-bond crystal 
(DVBC) states [19,14]. These DVBC states are charac- 
terized by the following nontrivial spin correlations: 

) ^ (5) 

when k approaches oo (for either all even i or all odd i 
only). Later, we present more evidence for the existence 
of DVBC Mott states; certain aspects of long wave length 
behaviors of DVBCs found for spin-one bosons are simi- 
lar to those of DVBC states in S* = 1 spin chains [20 - 25]. 

Valence-bond crystal states and more generally spin- 
Peirls states have been proposed in various models 
for strongly correlated electrons. Possible spin-Peirls 
states in 2D antiferromagnets (disordered) due to Berry's 
phases carried by hedgehogs were pointed out in [26] . In 
S't/fTVj-antiferromagnet models or their effective models, 
valence- bond crystal phases were discovered in various 
limits [27-29]. Most recently, valence-bond crystals in 
strongly-correlated systems due to Berry's phases in Ising 
gauge fields were addressed [30,31]. In quantum dimer 
models, crystal phases have been anticipated in [32] , and 
were reviewed in a recent paper [33] . 

Condensates for spin-one bosons with antiferromag- 
netic interactions in lattices also have fascinating prop- 
erties; most of arguments in the weakly interacting case 
can be carried out parallel to those developed for the 
continuous limit. Particularly, ground state wave func- 
tions should possess Ising gauge symmetries; and con- 
densates support interesting half-vortex-type topological 
excitations in additional to usual integer vortices [34]. 
Furthermore, the local spin dynamics is described by a 



nonlinear sigma model; for an individual condensate in a 
zero-mode approximation, a constrained 0(3)-quantum- 
rotor model can be introduced to characterize the spin 
dynamics. Due to spin-phase separation there are a vari- 
ety of novel spin correlated " condensates" ; particularly, 
we will argue that there should be spin singlet quantum 
condensates (SSQCs), in addition to usual polar conden- 
sates. 

Early indication for spin singlet quantum conden- 
sates, or rotationally invariant condensates came from 
the analysis of renormalization-group equations. The 
phenomenon of spin-phase separation implies that or- 
ders be established in either phase, or spin sector but 
not necessary in both. Particularly, there exist rotation- 
ally invariant superfluid states [34]. In fact, the failure 
of order parameters (when treated as quantum opera- 
tors) to commute with the Hamiltonian eventually leads 
to these fascinating many-body condensates of spin-one 
bosons. From this point of view, a polar condensate is a 
"classical" one, of which order parameter operators can 
be approximated as classical variables. 

The second relevant observation on spin singlet con- 
densates was made in [11]. There, the authors consider a 
limit where two-body repulsive interactions in ferromag- 
netic channels are much stronger than in singlet channels. 
In this limit, atoms form singlet pairs and condense; the 
resultant state is a condensate of spin singlet pairs. The 
solution in this limit is consistent with general arguments 
based on renormalization-group equations; furthermore, 
it indicates a microscopic realization of rotationally in- 
variant condensates suggested by renormalization-group 
equations. In [35], the authors further analyzed a few 
other spin singlet states, some of which break time reser- 
val symmetry and have fascinating edge properties. 

The most recent indication of SSQCs is the phe- 
nomenon of fractionalization in some Mott states. Con- 
sider each cold atom as a particle with spin 5=1 and 
"charge" Q = e. The spin-" charge" separation in Mott 
insulating states for spin-one bosons with antiferromag- 
netic interactions implies fractionalization of cold atoms 
in optical lattices; it also offers some insight into micro- 
scopic wavefunctions of fractionalized condensates. The 
ground state of interacting spin-one bosons in the strong 
coupling limit can be a DVBC; and this state supports 
either spinful but chargeless quasi-particles or charged 
but spinless quasi-particles in one-dimensional limit [19]. 
The low-energy sector of the Hilbert space can thus be di- 
vided into two subsectors of spin and charge. All the low 
lying excitations live in one of these two sectors. There- 
fore, each atom injected into a lattice naturally fraction- 
alizes into spinless charge-e {Q = 1, S = 0) excitations 
and chargeless but spinful {Q = 0,S = 1) excitations. 
Theoretically, the novel SSQCs appear naturally as a re- 
sult of condensation of spinless charge-e bosons. In this 
article, we substantialize those proposals made in [19]; 
particularly we demonstrate the microscopic structure of 
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possible chargc-c spin singlet condensates [36]. 

The fractionalization of cold atoms to be discussed 
in this article is an analogy of spin-charge separation 
proposed in strongly interacting electron systems. For 
strongly-correlated electrons in cuprates, the concept 
of spin-charge separation was emphasised in [37]; the 
phenomenon of spin-charge separation in doped anti- 
ferromagnets was studied afterwards in terms of effective 
compact gauge fields [38,39]. Many attempts have since 
then been made to understand the fractionalization of 
electrons in cuprates [27,32,40-46]. In some subsequen- 
tial works, fractionalization in quantum spin liquids was 
investigated in the context of effective Ising gauge theo- 
ries [28,29,47]. Most recent efforts on this subject can be 
found in [48,30,31,50,49,33,51]; especially, superconduct- 
ing states of "chargons" (charge-e objects) have been pro- 
posed in [30,51]. Some general discussions on spin-charge 
separation can be found in a remarkable book [52]. 

It is worth pointing out that from the point of view 
of fractionalized atoms, in some Mott insulating states 
only a fraction of each spin-one boson appears to be "lo- 
calized"; these states display rather interesting internal 
coherence as a result of " partial condensation" . For in- 
stance, nematic Mott insulating states can be effectively 
considered as condensates of fractionalized spinful but 
chargeless particles or spin-one spinous. In appendix C, 
we examine spin correlated Mott states from this point of 
view. Our discussions will be limited to cases of integer 
numbers of bosons per site; condensates for non-integer 
numbers of bosons per site will be addressed in a future 
paper. 

In this article, we develop a general approach for the 
studies of spin singlet states. First, we demonstrate that 
all SSMIs can be fully characterized by even- and odd- 
class quantum dimer models. In the second part of the 
article, we pursue the idea of SSQCs of spin-one bosons 
in lattices; particularly, we provide various evidence for 
these unconventional spin correlated condensates. In 
SSQCs, condensation takes place only in the charge sec- 
tor of the Hilbert space. SSQCs therefore exhibit dis- 
tinct long-range order. It is possible to distinguish them 
from pBECs in experiments because of remarkably dis- 
tinct macroscopic properties of SSQCs. 

In section II, we summarize results on spin- correlated 
Mott states of spin-one bosons in lattices and present mi- 
croscopic wavefunctions for some of the states. In section 
III, we investigate spin singlet Mott states for odd num- 
bers of atoms per site in one-dimensional lattices in the 
extremely small hopping limit. In section IV, we discuss 
excitations in spin singlet Mott states; we demonstrate 
the fractionalization of spin-one atoms in DVBCs in one- 
dimensional lattices. In section V, we present a qual- 
itative picture of SSQCs using a projected spin singlet 
Hilbert space; we argue that spin correlations strongly 
suppress one-particle hopping in an SSMI for an even 
number of bosons per site. In section VI, we introduce 



a fractionalized representation for interacting spin-one 
bosons in a lattice. In section VII, we briefly present re- 
sults on SSMIs and SSQCs in high dimensional lattices. 
In section VIII, we revisit low-dimensional SSMIs using 
an approach based on quantum dimer models; particu- 
larly, we generalize the results on one-dimensional lattices 
in section III to the entire Mott phase (one-dimensional). 
In section IX, we discuss SSQCs in one-dimensional lat- 
tices. In section X, we address the issue of the possible 
role of Berry's phases in one-dimensional lattices. 

It was demonstrated and emphasized in a few occa- 
sions that the problem of interacting spin-one bosons 
can be mapped into a constrained quantum rotor model 
(CQR) [34,11]. Technically, this allows one to acquire 
certain valuable intuition about many strongly-correlated 
states of spin-one bosons using the knowledge about 
the renormalization-group equations and topological field 
theories. The mapping is therefore a powerful device to 
study properties of spin-one bosons in optical lattices. 
The resultant CQR model will be our starting point for 
the studies of spin-one bosons in lattices. 

II. SPIN CORRELATED MOTT STATES: 
A SHORT REVIEW 

For spin-one bosons with antiferromagnetic interac- 
tions in (optical) lattices, the effective Hamiltonian can 
be conveniently expressed as 

Hlat. = -?E(V'LV'/a + h.C.) + ^ H + ^ - p,^ (6) 

{kl) k 

where the sum (kl) is carried out over neighboring sites 
k and I; in this section, we are only interested in high 
dimensional bipartite lattices. 

InEq.6, Eg = l/2/andi?c = 1/2C are two energy gaps 
in the excitation spectra of an individual condensate at 
each site, i is the hopping integral; /z is the chemical po- 
tential. Nq{'^ 1) will be the number of bosons per site. 
The dynamics is defined by the following commutation 
relations, 

[Sfc , ^If}] = ie°''^''Skltpk-f, [Pk, -^jj = ^klSafl'tpj^ (7) 

where 

= -ie«^^Vi^Vfe7> Pfe = 4a^ka (8) 

are the spin and "charge" operators for spin-one bosons, 
a(/3,7) = x,y,z. And 

[S^,Sf,]=ie«^^4,fe'S^,[S^,Pfe,] = 
[^kAi'04'0] = [^k,^k'0i^k'f}]=O (9) 

In addition. 
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Sfe = Pk{pk + 1) - tPka'^L'^kplpkp 



(10) 



where V'a^'a is a singlet pair creation operator. Following 

this identity, S| = Sk{Sk + 1) and Sk = 0,2,4, ...,2n, if 
Pk = 2n; otherwise Sk = 1, 3, 5, 2n + 1 if pfe = 2n + 1. 
For spin-one bosons with antiferromagnetic interactions, 
the ground state for an isolated site with pk = 2n atoms 
is a spin singlet and for = 2n + 1 is a spin triplet [5] . 

Therefore at low energies, the Hilbert space of the 
Hamiltonian is subject to a constraint that the sum of 
the number of particles pk and the total spin Sk, has to 
be an even number so that the many-boson wave function 
is symmetric under interchange of two spin-one particles. 
Say it differently, the Hilbert space displays the following 
selection rule 



(-1) 



Sk+Pk 



l(5'/c < Pk)- 



(11) 



This selection rule plays an extremely important role 
in both symmetry breaking states and spin singlet Mott 
insulating states [11,14,15,19,34]. 

In optical lattices, the hopping t can be varied contin- 
uously as a function of laser intensity; i decreases as the 
laser intensity increases. Es^c on the other hand depend 
on the scattering lengths ap, F = 0,2 



4irpo{a2 - Qq) 
3MNo 



.Er, = 



47rpo(2a2 + ap) 
ZMNo 



(12) 



We will focus on the case when Es is much less than E^. 
As laser intensities are varied, we then have at least three 
different situations; 



WcQR = -t 2^ rife • cos(xfe - xi) 

(ki) 



+ Y,EsSl + E. 



cPk - Pkl^- 



(15) 



And 



[Pk' , Xk] = ihSkk' ,Pk = ifi' 



9xk' 



(16) 



As the orientation of the 0(3)-rotor at site k is specified 
as rife, the total spin at each site is defined as the angular 
momentum of this rotor, following Eq.l5. The phase Xk 
further defines the orientation of an 0(2)-rotor; the num- 
ber operator corresponds to its angular momentum. The 
sum of Sk and pk again is subject to the constraint in 
Eq.ll. The Hamiltonian WcQR is locally invariant under 
a simultaneous inversion of rife and exp(ixfe). Introducing 
operators: 

6^^i...,Xk,---) = 'fi---,TT + Xk,---) 

rife,... )==*(..., -nfe,...) 
this invariance can be expressed as: 

{d^d^r'ncQnO^d^ = WcQR. (17) 

In the limit i), according to Eq.l5, a usual polar con- 
densate is established. A "mean field" polar condensate 
in this representation is 



*pBEc({nfe}, {Xfc}) oc Yl S{{nk ■ no)e 



i{Xk-Xo) _ 



1) (18) 



opt. 



» Er » Es 



n)E, » t{Vopt) » VEcEs-, 
m)Ec » s/EcEs » t{Vopt)- 



(13) 



In Eq.l3, t = Not is the effective hopping matrix element 
including a standard bosonic factor A^o- The wave func- 
tions of spin correlated states we are interested in, can be 
conveniently expressed in a coherent state representation 



(Xln. 



" V2!n ^ exp(-imx) 



m=No-5N 



|0), 



(14) 



with 6N <^No. 

So at each site, one introduces two " collective" coordi- 
nates rife and e''^'' , two unit vectors on a two-sphere S'^ 
and a unit circle S^ which characterize the orientation 
of 0(3) and 0{2) quantum rotors respectively. Following 
Eqs.6-11, the Hamiltonian Hut. {No, SN » 1) in this rep- 
resentation can be expressed as Hcqr of a constrained 
quantum rotor (CQR) model 



where all directors rife point in the direction of no and 
all Xk are "locked" at the same value of xo- The wave 
function in Eq.l8 explicitly exhibits the local Ising gauge 
symmetry emphasised in [34], i.e. 

OfeO^fpBEC = fpBEC. (19) 

Microscopically, 

represents a condensate at a state with zero crys- 
tal momentum Q = 0, along axis Hq; Xga ~ 
-y= X]exp(ifc(5)Xfe„, here Vt is the number of sites in 
a lattice. 

In the limits of ii) and iii), particles are localized at 
each site and there are Nq atoms per site. The wave- 
function can be factorized into 

^'({nfe}, {xfe}) = ^NoiWk}) n exp(iiVoXfc) 

^ v27r 



(21) 
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FIG. 1. Ordering of two order parameters, unit vectors n 
(a light arrow defined on a two sphere S^) and exp(ix) (a dark 
arrow defined on a unit circle , drawn in the YZ-p\ane of 
the two sphere) in different correlated states. a)pBEC, b)NMI 
and c) SSMI. Two spheres and unit circles (5*^ and S^) for a 
few different lattice sites are shown here. 

In terms of collective coordinates n^, the effective 
Hamiltonian for Mott insulating states can be expressed 
as 

Ti-Mi = —Jex ^ Qa/3(nfc)QQ/3(n;) + Eg ^ S^; 

(kl) k 
Qapin) = UaUp ~ ^Sal3. (22) 

The exchange interaction J^x is of the order of t^Ec^^. 
The Hilbert space observes the following symmetry under 
a local inversion O", 

0,"*A.„ = i-if^^No; (23) 

Finally, the Hamiltonian is invariant under the local in- 
version 

(6^)-1Hmi6" = Hmi- (24) 

In the limit of ii) , one can identify the ground state as 
a nematic Mott insulator. In this intermediate regime, 
the director n is localized on a two-sphere S^. The wave- 
functions for an even (e) and odd (o) number of atoms 
per site are 

nMii{nk})^YlS{in,-nor-l) 

k 

*NMi({nfc}) « n '^(("'^ • no)' - • "o)- (25) 

fc 

One notices that 



so that the constraint in Eq.23 is satisfied. Eq.25 indi- 
cates that at each site all localized atoms condense in an 
identical spin-one state characterized by i/'lnoa; micro- 
scopically (see Fig. 2a)) 

*-j«n^^i^^i^io) (27) 

which corresponds to a maximally ordered state. 

A direct calculation of order parameters defined in 
Eq.l in section I yields 

pBEC -.Oi = y^Von^ exp(zxo); 
NMl:Oi^0,Olp^ No{n^nfj~ha,f3). (28) 

Nematic Mott insulating states can be observed at tem- 
peratures lower than the exchange energy Je^.. 

Nematic states break the rotational symmetry. The 
internal space and its first homotopy group are 

5-2 

7^ = — ,7^l(7^) = ^2. 

^2 

These states support interesting 7r-spin disclinations 
(some general discussions on the effects of Ising symme- 
try on atomic defects can be found in an early work [34] ) . 

In the limit of iii), n and x s-re delocalized on the unit 
sphere and unit circle respectively; for an even 
number (Nq = 2n) of particles per site, the wave function 
for a spin singlet Mott insulator (SSMI) can be written 
as 

*isMi({nfc})«n^o^o("fc)' (29) 
fc 

Yo.o is the zeroth spherical Harmonics. ^'gsMi rota- 
tionally invariant unlike NMI or pBEC; the microscopic 
wave function which yields the desired collective behavior 
specified in this limit is (see Fig. 2b)) 

^^^^ Y V(2"+l)- 

The situation for an odd number of atoms per site 
is more involved and will be discussed in the next sec- 
tion. As mentioned, more detailed discussions on quan- 
tum spin ordered or disordered phases can be found in 
two very recent works [14,15]. 
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(b) 
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FIG. 2. Schematic of microscopic wave functions of an NMI 
(a), an SSMI for iVo = 2n (b) and a DVBC for A^o = 2n + 1 
(c). A dot carrying a double-headed director pointing in di- 
rection n represents a spin-one boson in a state ^^nQ,|wac); 
each pair of dark and light dots connected with a ring stands 
for a spin singlet pair of two spin-one atoms. 



III. DVBCS IN ONE-DIMENSIONAL MOTT 
INSULATING STATES: SMALL HOPPING LIMIT 

Consider one-dimensional Mott states when the hop- 
ping is much less than Ec^s- For an odd number (A'o = 
271 + 1) of atoms per site, each two unpaired atoms at 
two neighboring sites form a spin singlet and the ground 
state is a dimerized- valence-bond crystal state [19]. The 
previous results suggest the following wave functions (see 
Fig. 2c)) 



*SSMI 



n 



xjc/3X(fe+i)^|0). 

(31) 



(2n + 3){2n+l) 



The product Y[' is carried over all odd sites (fc — 2m + 1) 
or all even ones (fc = 2m). Since the key idea of SSQCs is 
intimately connected with properties of DVBCs, in this 
section we will present evidence for DVBC states in one- 
dimensional lattices. 

Though dimerization occurs in all one-dimensional 
Mott states for an odd number of particles per site, we 
start with a simple limit when 



t <C V EcEs, Ec] 



(32) 



we would like to demonstrate the origin of these dimer- 
ized states of spin-one atoms from the point of view of 
bilinear-biquadratic (BLBQ) spin models. Discussions 
on more general cases can be found in section VIII. The 
Hilbert space at each site is defined in Eqs.7-11. For 
odd numbers of bosons per site(Afo S> 1), the low energy 
Hilbert space for each isolated site is spanned by states of 



total spin S = 1,3,5.... Therefore, as Es is much larger 
than one can further truncate the Hilbert space at each 
site into a space of a spin-one particle (see FIG. 3). 

The effective exchange interactions between two adja- 
cent condensates can then be obtained by examining the 
excitation spectrum of two coupled sites. At t = 0, the 
ground state of two decoupled condensates has a nine- 
fold degeneracy. At any finite t much less than Es^c: 
because of virtual hopping between two sites, the nine- 
fold degeneracy in the truncated space is lifted and the 
resultant spectrum is 



E[S) = -asJ, 5== 0,1,2; 
as calculated for two coupled condensates are 
71 (5) , 12{S) , 73(5) 



as 



and 



1 - 4a 
1 



1 + 2cs 



1 8c, ' 



(33) 
(34) 



g,73(0) = |,73(l) = ^; 
2 7 



71 (0) 
72(2) 

72(0) =7i(l)= 72(1) -7i(2) = 0. 



72(2) = —,73(2) = 



(35) 



Cs — Eg/Ec] and ao is always larger than a 1,2. 

The truncation can be carried out in a similar way in 
a lattice. The following reduction in dimensions of the 
Hilbert spaces takes place when the hopping is weak 



,(iVo + l)(A^o + 2). 



Vt 



Vt = 3 



(36) 



Here Vr is again the number of sites in a lattice and the 
number of bosons at each side No is much larger than 
one. In the truncated Hilbert space of dimension Vt, 
Eq.33 suggests the following effective Hamiltonian 



{kl) S 



(37) 



where VsiSk, S/) is a projection operator for a total spin 
S state of two condensates at neighboring sites k and I 
(S| = Sf = 2). 



H 



MI 



H 



bb 



B 



.S=l 



B 



FIG. 3. The truncation of the Hilbert space of two conden- 
sates A and B (A^o = 2n -I- 1) when t is much less than E^.s- 
The effective Hamiltonian in the truncated space is Tib.b.. 
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For onc-particlc per site (A'o = 1) at small hopping 
limit, a straightforward calculation yields values of as, 
5 = 0,1,2, 



ai = 0, ao 



l-2c, 



■,a2 



1 + Cs 



and again ao > a2 > oi- 

Now we can proceed to one-dimensional lattices where 
the situation is relatively well understood. For an odd 
number of bosons per site {Nq = 2k + 1), the Mott insu- 
lating state should be a DVBC with a twofold degeneracy 
[19]. To obtain this result, we take into account the pe- 
culiar property of the truncated Hilbcrt space and energy 
spectrum. Following Eq.37, one shows that the problem 
of interacting spin-one bosons with an odd number of 
bosons per site can be mapped into the BLBQ spin-one 
spin chain model. 



'^b.b. 



= cosTj^^ Sj • Sj + sin 77 ^^(S, • Sj)^ 



(38) 



if T] satisfies tan 77 = (ai — ia2 — |q:o)(q!i — 0:2) ^ and the 



sign of sin 77 is chosen to be the same as ai — — fao- 
Cg varies from to 1/4 when iVo = 2fc + 1 ^ 1, and from 
to 1/2 when Nq = 1; consequently, r] varies as 



-->??> ?7o. 



(39) 



The value of rjQ in general depends on models used in dif- 
ferent limits; for the large- iVo case { Nq = 2k + 1 ^ 1), 
r/o = — f — arctan i while for A^o = f, Vo = — 37r/4. For 
parameters given in Eq.39, following discussions in [40] 
we arrive at conclusions that the ground states should be 
dimerized valence bond crystals or DVBCs. The wave 
function of DVBCs for A^o = 2A; -|- 1 atoms per site is 
proposed at the beginning of this section. At 77 = 0, the 
model describes a Heisenburg antiferromagnet; solutions 
in this limit were also studied in [53]. 

Taking into account Eqs.8,31, one obtains spin-spin 
correlations in DVBC states. The DVBC state breaks 
the crystal translational symmetry and has a twofold de- 
generacy; one of degenerate states is characterized by the 
following correlation functions (m ^ n) 

{Sm ■ S„) CX -i(l - (-1)™) (1 + {-ir)Sn,m+U 

{Sm ■ S™+lS„ • S„+i) OC ^(1 - (-1)'")(1 - (-1)"). 

(40) 

Because of the selection rule discussed in Eq.ll in 
section II, for even numbers of bosons per site, the 
Hilbert space is spanned by states with total spins of 
S = 0,2,4..., Nq. The lowest energy space can be trun- 
cated into a spin singlet state. As the coupling t between 
sites is much smaller than Eg, the exchange interaction 



is an irrelevant operator; the groimd state remains to be 
a non-degenerate spin singlet Mott insulator as far as t 
is much smaller than E^. In this limit, spin-spin correla- 
tions fov n are zero, i.e. 

(S„ • S„) « 0. (41) 



IV. EXCITATIONS IN SSMIS 

In this section, we will demonstrate; charge-e and 
charge-2e excitations. The existence of spinless but 
charge e (5 = 0, Q = 1) excitations which do not carry 
full identities of spin one-bosons implies the fractional- 
ization of spin-one atoms. These spinless objects also 
provide important hints on microscopic wave functions 
of SSQCs discussed in the next section. 



A. Charge-e excitations in one-dimensional DVBCs 

In terms of creation operators of a valence bond at link 
77 which connects two neighboring sites i and j 

<l>\ = ^(hl,ih],-i + hl.-ihii - hl,oh],o), (42) 

the wave function of one-dimensional DVBC states can 
be written as 



*DVBC = 11'^^ 



(43) 



The product H' is carried over all even or odd links; ht ^ 
is a creation operator of a condensate at site i with total 
spin S = 1, Sz = m. 

A one-dimensional DVBC state supports two kinds of 
elementary excitations which are of topological nature 
[19]. A spinless but "charge-e" excitation (Q = 1,5 = 0) 
represents a spin singlet state with an extra atom com- 
pared with the Mott ground state; a chargeless but spin- 
ful excitation (Q = 0,5 = 1) on the other hand is a 
spin-one state but with the same number of atoms as the 
Mott state. Both excitations are kink-like and created 
only by nonlocal operators. 

An S = l,Sz =171 excitation can be created by the 
following product operator 



7,m 



V'r. 



G.G.^\,m 



(44) 



rieCy 



The product is carried over all links 77 along path Cj 
which starts at site 7 and ends at infinity; furthermore, 
Cj is chosen such that the first link along the path is occu- 
pied by a valence bond. Pq ^ is a generalized Gutzwiller 
projection to project out states with 71-particles per site. 
Therefore, Eq.44 indeed represents a spin-one domain 
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wall soliton {m{Sz) = 0,±1) located at site 7 in the 
DVBC state as shown in Fig. 4a). 

The hopping matrix of domain-wall along an one- 
dimensional lattice is 



+2J 



Taking it into account, we obtain the following band 
structure for spin-one excitations; 



E{Qx) = Jao(l - -cos2Q^) 



(45) 



with Qx defined as a crystal momentum of excitations, 
— 7r/2 < Qx < tt/2. These excitations are purely mag- 
netic and involve no extra bosons; therefore are "neutral" 
(Q = 0). 

Besides spin-one "neutral" excitations, there are also 
"charged" spinless {S = Q,Q = 1) domain wall excita- 
tions. A "charged" excitation with a positive or negative 
charge {Q — ±1) is defined by a creation operator i?^ _|_ 
and 



7,m 7,- 



(46) 



fpi mi''l'i,m) is a creation (annihilation) operator for a 
spin-one atom with Sz = m. The situation in one- 
dimensional optical lattices is therefore similar to con- 
ducting polymers [54]. 

The splitting of the low energy Hilbert space suggests 
the fractionalization of atoms. It also suggests a possible 
quantum condensate of fractionalized particles of charge- 
e. Following Eq.46, as it is added to an one-dimensional 
DVBC, a spin-one atom fractionalizes into a spin-zero 
domain wall (with an extra atom thus " charged" ) and a 
spin- one "neutral" domain wall because of the twofold 
degeneracy for A^o = 2n + 1 atoms per site. 



(a) 



(b) 




FIG. 4. Microscopic wave functions of excitations in 
one-dimensional optical lattices, a) A spin-one kink exci- 
tation (S = 1, Q = 0) in a DVBC; b) A charge-e excita- 
tion (Q = 1, 5 = 0) in a DVBC; c) A charge-2e excitation 
(Q = 2, 5* = 0) in an SSMI for an even number of atoms 
per site. In d), we also show an S = 0, Q = 1 state in an 
SSMI {No = 2n); creation of this state results in a string of 
singlet bonds and is an energetic catastrophe. Locations of 
excitations are indicated by dashed lines. 



B. Charge-2e excitations in one-dimensional SSMIs 

The only spinless excitations in one-dimensional SS- 
MIs for an even number of atoms per site are charge-2e 
ones. To demonstrate it microscopically, one considers a 
limit when Eg is infinity. All singlet pairs are bound and 
elementary charge excitations are charge-2e objects. At 
large but finite Eg, one finds at each site, the only spin- 
less excitations are created by adding or removing singlet 
pairs with Q = 2 and 5 = 0. These again correspond to 
charge-2e excitations. 



A PROJECTED SPIN SINGLET HILBERT 
SPACE BASED APPROACH 



To apprehend the key physics in SSMIs and SSQCs, 
we first provide an intuitive approach to this problem. In 
a projected Hilbert space, we are going to demonstrate 
that the problem of spin-one bosons with antiferromag- 
netic interactions can be mapped into spinless bosons 
interacting with Ising gauge fields. 

Let us consider a projected spin singlet Hilbert space 
which satisfies the following conditions: 

a) each atom has to form a spin singlet with another 
one either at the same site or a nearest neighboring site; 

b) each link has to be occupied by at most one spin 
singlet pair of atoms (a valence bond). 

c) Following a), b), the parity of the number of valence 
bonds emitted from a lattice site is the same as that of 
the number of particles at that site. For instance, for a 
site with one particle, one of the links connected to that 
site has to be occupied by a valence-bond. For a site with 
two particles, either two or none of the links are occupied 
by valence-bonds. 

In Fig. 5, we show examples of states in this projected 
Hilbert space. We employ a valence bond counting oper- 
ator 



1 



ife,fc+i — 



'k,k+l 



(47) 



at each link rj = {k,k + 1), the eigen values of which are 
dk,k+i = 0, 1. So we define the projected spin spin singlet 
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Hilbert space as the one constructed out of eigenstates 
of {pk} and {dk,k+i} 



\...Pk,dkM+uPk+i-), if (-ir+s.^^" = 1. 



(48) 



Here is the number of particles at site k, dk,k+i is 
the number of valence bonds at a link connecting site k 
and fc + 1. The local constraint in Eq.48 follows point c) 
discussed above; rj represents a link connected to site k. 





(e) 



(f) 




(g) 



FIG. 5. Examples of states in a projected spin singlet 
Hilbert space, a), b) and d), e) are states in this space. States 
in c) or f) and g) can be "deformed", respectively, into states 
in a) or d) and e) by locally breaking valence bonds and con- 
structing intra-site singlets without involving atoms at a third 
site. States in c), f) and g) have higher energies than their 
counterparts in a), d) and e) and do not belong to the pro- 
jected Hilbert space. In a)-g) dots connected by rings are spin 
singlet pairs of spin-one bosons. 

Consider hopping of an atom between two neighbor- 
ing sites fc, I (an even number of atoms per site in initial 
states, see Fig. 6a). 



XL 



f (^fcai^/cai) 



(49) 



v/(2«-l)! , 

In a projected Hilbert space, hopping of atoms sim- 
ply leads to transitions between the following states 

\---Pk, dk,k+l, Pk+1---) 

\---Pk,dk^k+l, Pk+1---) 

{■■■Pk - 1, 1 T dk^k+i,Pk+i 



1.. 



(50) 



for dk,k+i = 1 or in initial states. One obtains an ef- 
fective description of hopping of spin-one atoms in the 
projected Hilbert space. 




m 1^ 



+1 



(a) 



+1 



(b) 



FIG. 6. Hopping of atoms in spin singlet Mott states and 
its effective description. For an even number of particles per 
site in a), hopping of an atom adds a singlet bond between 
two sites; in the lower part of a), this process is represented 
by changes in "charges", at each site (numbers below 

dots) and "countings" of valence bonds d,, at each link (num- 
bers above links), b) is for an odd number of particles per 
site. Notice that the parity of -I- foj^bfc at each site is 

conserved when atoms hop around; 77 is a link connected with 
site k. The description in a) and b) leads to an effective the- 
ory on spin-one atoms in ID lattices. In a) and b), we only 
keep states in the projected Hilbert space introduced in this 
subsection. 



Moreover, at zero hopping, energies of states with dif- 
ferent particles are given by "charging" energies in Eq.6. 
And to create valence bonds one needs to break intra site 
singlet pairs; and thus the energy cost of a valence bond 
is Es if Jex is much smaller than Eg. Therefore in the 
projected Hilbert space \...pk,dk,k+iT Pk+i^^^)-, we obtain 
the following mapping 



Pi 



^ f,t f, Pk 

^ 2T,4 

21 



(51) 



(using as a linear combination of raising and lower- 
ing operators defined with respect to d). The resultant 
Hamiltonian is the same one as in Eq.103, which will be 
derived in a more formal way. It is easy to confirm that 

Tb ~ Es when Eg > Jea:- 
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t ° • • ^ > 

2n+l 2n+l ^^-^ 2n+l 2n+l 2n+l 

2n+l 2n+2 ^-^^ 2n+l 2n+l 2n+l 

^ ^ ^ .1 ^ 1 ^ 

2n+l 2n (c) 2n+l 2n+2 2n+l 

2n 2n 2n 2n 2n 
(d) 

• ° ^ ^1 ^ ^ -1 ^ 

2n 2n+l 2n 2n 2n 
(e) 

> ° • • m ' m 

2n 2n+l 2n 2n-l 2n 
(f) 



This difference between lattices witfi even or odd num- 
bers of particles per site mentioned above plays a remark- 
ably important role when hopping is taken into account. 
In Fig. 7, we illustrate hopping of atoms schematically 
from this point of view. Hopping of atoms in an " odd" 
lattice results in kink-anti kink pairs which are weakly 
interacting. However, in an "even" lattice hopping is 
strongly suppressed because of strings of valence bonds 
created between particle and hole excitations as shown in 
Fig.7 f); it implies that only pairs of atoms or charge-2e 
objects propagate along the lattice in this case. In the 
following sections, we implement this idea more quanti- 
tatively using some properties of Ising gauge fields. 

Similar analysis about one-particle hopping in spin sin- 
glet states can be easily carried out in high dimensional 
lattices. For even numbers of atoms per site, one arrives 
at the same conclusions as in one dimensional lattices. 
For odd numbers of atoms per site, using the mapping 
argued in this section, one can also investigate hopping 
in possible DVBC states in square lattices (see Appendix 
A). 



FIG. 7. Schematic of one dimensional SSMIs and hopping 
of bosons in these states, a) a DVBC for an odd number 
of particles per site; b) an 5* = and Q = 1 kink excita- 
tion in a DVBC; c) hopping of an atom results in a kink-anti 
kink pair; d) an SSMIe for an even number of particles per 
site; e) an 5* = and Q = 1 "excitation" in an SSMI (for 
N — 2n) emitting a string of valence bonds; f) hopping of 
a boson is suppressed in an SSMI (for A'^ = 2n) because of 
a string of valence bonds between particle-hole excitations. 
Numbers appearing at each link are countings of d and those 
below dots are numbers of particles at each site; light dots 
are "charged" and along light links valence bonds have been 
removed or added to ground states. 

The effective Hamiltonian is convenient for the study 
of spin singlet states. Let us consider spin singlet Mott 
insulating states discussed in section II, III; In the pro- 
jected space, when Es ^ Jex, ground states can be read 
out as 

Pk = 2n + 1, 4,fc+i - ^(1 ± (-1)'=), for No = 2n + 1; 
Pk = 2n, dk^k+i = 0, for A^o = 2n (52) 

as illustrated in Fig. 7a) and 7d). 

The effective Hamiltonian in the projected spin singlet 
Hilbert space is also particularly useful for the discussion 
of hopping of atoms across lattices. Following the con- 
straint in Eq.47, the insertion of an atom to site k changes 
Pk by one unit and leads to an = 0, Q = 1 kink-like 
excitation (see Fig. 7b)); this reproduces an spin singlet 
charge-e excitation in Fig. 4b). On the other hand, for an 
even number of particles per site, adding an atom to a 
site results in an S" = and Q = I string-like excitation 
(see Fig.7 e)), this was also discussed in Fig.4d). 



VI. A FRACTIONALIZED REPRESENTATION 

For a discussion on SSQCs or more general fraction- 
alized quantum condensates (FQCs), we find that it is 
most convenient to introduce a representation involving 
"chargeless" {Q — 0) spin-one operators and "spin- 
less" {S = 0) charge-e(Q = 1) operators 6; 

4a = hxL- (53) 
a, b, p and S satisfy the following algebras 
[Pk,b\] = ^kibl, [pfe,a|] = 

[S^,a|^]=»e"^^5fe;al^,[S^,6|]=0. (54) 

Eq.53 suggests that 

Pk = b\bk, SfcQ = — ieQ^^aJ.^afe^; 

the spin density is carried by a-particles, or spin-one 
spinons and charge density is carried by 6-particles, or 
chargons. It implies a possibility of spin-" charge" sepa- 
ration in cold atoms. General discussions on spin-charge 
separation in strongly correlated electrons can be found 
in [27,28,30-33,37-49]. 

Obviously, the enlarged Hilbert space spanned by {a^} 
and {b^} includes unphysical states violating the sym- 
metry of many-boson wave functions. Because of the 
symmetry constraint, all physical low energy states have 
to display a superselection rule discussed in details in 
section II. Namely, for an even number of atoms A"o, 
the low energy Hilbert space is occupied by states with 
5 = 0, 2, 4, 6..., Nq while for an odd A'o: the Hilbert space 
is spanned by states with S = l,3,5...,Aci. 
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Following the identity in Eq.lO, if one has an even num- 
ber of a'^-particles, resultant total spins are S = 0, 2, 4... 
and for an odd number, spin states are those with S = 
1, 3, 5.... To produce desired symmetric states of spin-one 
atoms, one has to require at least that the total numbers 
of a- and b- particles, or spinous and chargons, be even. 
If the average numbers of a- and b- particles at each site 
are much larger than one, it is possible to generate cor- 
rect low energy Hilbert spaces by imposing the following 
constraint so that for every state 

Cfe* = ^f, Cfe = exp ( ± inial ■ afe + 6^6fe]) • (55) 

The constraint in Eq.55 is precisely to exclude the un- 

physical states in the enlarged space spanned by and 
ajj which violate the selection rule in Eq.ll. 

Considerations based on previous works on discrete 
gauge symmetries [34,11] indicate that the Hamiltonian 
in Eq.6 should be equivalent to the following one in this 
fractionalized representation, 

= E If + ^ - ^fc^ - *E i^a^i-bibi + h-c), 

k (kl) 

(56) 

In appendix B, we provide further evidence for the equiv- 
alence between Eq.56 and Eq.6. 

Let us emphasize that we are interested in the follow- 
ing limit where the equivalence has been established; 

{pkb) = No, (pka.) = iVa (57) 

and both A^o and A^a(^ ^o) are much larger than unity. 

Pkb = blbk, Pka = ^ia^ka, 

No is the number of spin-one particles at each site deter- 
mined by the chemical potential; the exact value of A^a 
depends on the trucation of the low energy Hilbert space 
and is not important for the rest of discussions. 

The resultant Hamiltonian is invariant under a local 
gauge transformation 

C^^HpRCk = Ti-FR- (58) 

We should also emphasis here that only bilinear oper- 
ators of a^ , are invariant under the local gauge trans- 
formation defined in Eqs.55,58, 

^k ^^ka^kffik = ^ka^kp^ 

^k ^^L^I^fe = ^L'^L 

Ck%blCu=blbl 

(59) 

On the other hand, all linear operators transform non- 
trivially under this transformation; 



^k ^^ka^k — -a|,^, 

Ck%Ck = <■ 

(60) 

From a standard point of view, these linear operators 
carry charges defined with respect to the gauge trans- 
formation in Eq.55 while bilinear operators in Eq.58 are 
charge neutral. 

This property of "fractionalized particles" is more ex- 
plicit if we examine the corresponding action of the 
Hamiltonian derived in appendix B. Besides spinous and 
chargons, we also have introduced discrete gauge fields 
defined at links of the lattices [55-58] . The constraint 
discussed above results in Ising gauge fields a^ki = ±1 in 
D + 1 Euclidean space. And 

^ki = Y^i^ki-'^ki)^ 

■^^-^kU^k'l'] = ^Skl,k'l'(^ti- 

Following appendix B, an inversion of n ^ — n, or 
exp(ix) — cxj>(ix) corresponds to a discrete gauge 
transformation a —a, or 6 ^ —b. The action is there- 
fore manifestly invariant under the following local Ising 
gauge transformation 

bl ^ Ctkbi, 4 ^ Ofea^, a'ki ^ Ofea^; (61) 

with Ofc = ±1. 

As indicated here, a^ and b^ are indeed matter fields 
carrying charges of Ising gauge fields. Taking into ac- 
count minimal coupling between gauge fields and matter 
fields, alternatively we construct a Hamiltonian which 
explicitly involves gauge fields ct|(, spinon fields a^ and 
chargon fields b]^. 

g2 2 

k 

""kl (aL^ia + blbi + h.c) , 

{kl) 

(62) 

The Hilbert space of the Hamiltonian in Eq.62 is again 
subject to the constraint in Eq.55. The Hamiltonian is 
invariant only under the following generalized transfor- 
mation 

C'^=4®exp(±i7r^l^), 
+ 

{C'^)~^HmgCI, = Hmg- (63) 

The sum over the cross (-I-) is carried over all links con- 
nected with site k. Calculations indeed suggest that the 
low energy physics of the new Hamiltonian be the same 
as the one in Eq.6 [59]. Both Eq.56 and Eq.62 will be em- 
ployed for discussions on FQCs and specifically SSQCs. 
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VII. SPIN SINGLET STATES IN HIGH 
DIMENSIONAL LATTICES 

To facilitate discussions on SSMIs and SSQCs or more 
general FQCs, we introduce the following correlation 
functions 



g\k,i) = {blbi)^ 



(64) 



evaluated at a fixed gauge similar to a "minimal gauge" 
suggested in [57]. In this gauge, the number of links 
where cr|; = —1 is minimal for a given distribution 
of gauge fields. The two-point correlation functions in 
Eq.64 are expected to diagnose possible condensation in 
the spin or charge sector of the Hilbert space. In addi- 
tion, we introduce local pairing order parameters 



A\k,k) = {bibi) 



(65) 



which are manifestly gauge invariant under the local 

gauge transformation defined in Eq.58. 

We further assume that the two parameters 

IP TP 

-tils Hjc 

can be varied independently. Two critical values Gsc, Gcc 
are introduced below to facilitate discussions on SSQCs 
and SSMIs. For the purpose of demonstration, in this 
section we limit ourselves to lattices with an even num- 
ber of bosons per sitc(iVo = 2n), though a generalization 
using the scheme developed in section V is possible. 

Large Hopping limit: Gs > Gsc, Gc > Gcc Under 
this condition, both spinous and chargons condense and 
the ground state represents a polar condensate or pBEC. 
The gauge fields are screened and cr^ — 1 at each 
plaquette; here □ represents an elementary plaquette. 
There are two branches massless spin wave excitations in 
this rotational- symmetry-breaking state and one branch 
massless density mode. The ac Josephson oscillation fre- 
quency is precisely the difference between the chemical 
potentials of two condensates. 

Intermediate Hopping limit I: Gs > Gsc, Gc < Gcc 
When Es is much less than Ec , this corresponds to a limit 
where EcGcc > t > EsGsc- Following the Hamiltonian, 
chargons are completely depleted from the condensate. 
We are particularly interested in a self-consistent solu- 
tion where spin-one spinous "condense" one way or the 
other. States in this limits are incompressible and there- 
fore naturally correspond to Mott insulating states; con- 
densation of spinous leads to nematic long range order. 
Nematic Mott states were studied briefly in section II; 
in appendix C, we provide an alternative description on 
NMIs using this fractionalized representation. 



Intermediate Hopping limit II: Gs < Gsc, Gc > 
Gcc This limit can be relevant to atomic gases in high di- 
mensional optical lattices, especially when Eg is not too 
small compared with Ec- We are going to present results 
in this case, for the completeness of analysises on SSQCs 
and more important to facilitate discussions on low di- 
mensional optical lattices in the next section. This case 
represents a situation where EsGsc > t > EcGcc- 

Discussions in this case can be carried out parallel to 
the previous case analyzed in appendix C. All spinful 
excitations are fully gapped; after integrating out these 
excitations, we obtain an effective Hamiltonian in terms 
of chargons 



'2 

k {kl} 



(kl) a a 



(66) 



Fb, are again functions of t, Es- Ksb is much less than 
Fb when Gg is small, 



^ « 1(F6 ~ &Es); 

I 6 



(67) 



and expected to be divergent close to Gsc- The induced 
constraint on the Hilbert space is 



C-f * = = exp {m[blbk ' ^ ^ 



-I- 



^]). (68) 



And the Hamiltonian is also locally invariant under the 
action of C-, 



fb 

k ' 



(69) 



If Ff, is much less than Ksb, chargons condense and the 
phase coherence is maintained, 



g\k, A; + oo) ^ 0, gip{k, k + oo)=Q 



(70) 



and the rotational symmetry is luibroken. This state is 
compressible but should have zero spin susceptibility at a 
low temperature limit. We define this state as a charge-e 
SSQC. 

One can also carry out discussions on a limit when t is 
much less Eg, or 



K, 



» 1. 



sb 



(71) 



And gauge fields are expected to be unity at each link 
over the whole lattice; under the constraint in Eq.55, the 
reduced Hilbert space is subject to a constraint 

Cl^ ^^,Cl=eM^^blbk). (72) 

Thus, the low energy physics is determined by the pair 
hopping of chargons; 
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k kl 

{Cly^TipbCl = Hpb- (73) 

and 

t2- ^{Tb'-^eE,). (74) 

i b 

When t2 = N^ii ^ E^, the ground state is a condensate 
of paired chargons; in terms of chargon-operator b\ 

I52) =n^^oxy.cxp(<i>J(0))|0), 

k 

$t(Qo) = Ao5]6tj^+/_,. (75) 

The projection operator VnoxVt acting on a general- 
ized BCS pairing wave function is to project out NqVt- 
particle states; A'o and Na arc even numbers. Again, 
Gaff{k, I) and G''{k, I) vanish as — / approaches infinity 
but Eq.75 indicates 

A''(fc,fc) = Y#^ (76) 

in the ground state. Finally, wc should emphasize that 
this charge-2e SSQC could be alternatively considered as 
condensation of " charged" soliton-anti soliton pairs from 
a point of view of section IV. 

Small hopping limit: Gg < Gsci Gc < Gcc 

We are going to revisit the strong-coupling limit of 
spin-one atoms in Mott insulating states: wc will em- 
ploy the fractionalized representation to express the wave 
functions of SSMI. Discussions are limited to even num- 
bers of bosons per site. 

In a strong-coupling limit Gs,Gc "C 1 and ti <C Es, 
t2 <C Ec, gauge fields of ct^ "condense". The Hilbert 
space of spinons and chargons explicitly satisfies a con- 
straint C^C^* = (C^ is defined in Appendix C.) 

As t is much smaller than the interaction energies, the 
Hamiltonian commutes with local spin and number op- 
erators in the leading order of Ec^s', 

[n,sl] = [n,pl]^o. (77) 

Therefore, the ground state should be an cigcnstatc of 
these two local operators and the wave function can be 
expressed as a projected BCS state of chargons 

l53)=P^G.l52) (78) 

where Vg^q. ^ generalized Gutzwiller projection for 
2n-cliargons per site. One can show that S*^ = for 
any k and Eq.78 represents a spin singlet Mott insulator 
(SSMI) studied in section II. 

Before turning to SSQCs in one-dimensional lattices, 
we first revisit SSMIs in one-dimensional lattices to fur- 
ther understand the even-odd effect discussed in section 



III. Wc are going to generalize previous results on one- 
dimensional SSMIs using mappings to an instanton gas 
model of Ising gauge fields and quantum dimer models. 



VIII. SSMIS IN LOW DIMENSIONAL LATTICES: 
A GENERAL APPROACH 

A. Instantons in Ising gauge fields 

From the point of view of Ising gauge fields introduced 
in the previous section, distinct properties of Mott states 
for odd and even numbers of bosons per site discussed in 
section III can also be attributed to a topological term in 
the effective actions for even and odd numbers of bosons 
per site. In this section, wc are going to revisit SSMIs in 
low dimensional lattices employing the effective descrip- 
tions introduced in section VI. We will reproduce results 
in section II and III; furthermore, in one-dimensional lat- 
tices we are going to generalize; the previous results on 
DVBC states in the limit t <C y/EgEc (see section III) to 
the entire Mott phase. 

The quantum problem of the Hamiltonian Hfr, or 
equivalently the Hamiltonian of constrained quantum ro- 
tor model HcQR defined in d dimensions can be mapped 
into the following model in (d+l) dimensions (appendix 
B), 

nFR< ^WcQii^ 0(2)® 0(3)0^2. (79) 

The 0(2) (g) 0(3) non-linear sigma model characterizes 
the dynamics of two variables n, exp(zx) introduced in 
section II, or spinons a+ and chargons fe+ in discussed in 
section VI. 

In high dimensional lattices when t is much less than 
Eg and Ec, spinons don't condense and chargons are lo- 
calized because of repulsive interactions; and the direc- 
tor n becomes disordered as suggested by a renormaliza- 
tion group equation analysis on the 0(3) nonlinear sigma 
model. For an even number of atoms per site, this strong 
c;oupling fixed point can be confirmed by directly solving 
the Hamiltonian Hmi in the same limit in all dimensions. 
In one-dimensional lattices, spin-one spinons again do not 
condense and the director n is quantum disordered at any 
finite Eg because of long range fluctuations. 

For an odd number of atoms per site in the absence of 
hopping, the low energy space is highly degenerate (see 
section HI) and the situation is more delicated. On the 
other hand, from the point of view of effective actions, the 
difference in numbers of atoms per site is only reflected 
in a Berry's phase term(see Appendix B). As far as low 
dimensional lattices are concerned, we speculate that this 
topological term doesn't alter the renormalization group 
flow of the 0(3)-model and arrive at the conclusion that 
for an odd number of atoms per site, spinons also do 
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not condense and the director n is also disordered as the 
strong-coupUng hmit is approached [60]. This observa- 
tion at least can be justified in one-dimensional lattices 
by comparing the results obtained below and the ones 
presented in section III. 

In this paper, we would like to assume a quantum dis- 
ordered phase in square lattices for odd numbers of atoms 
per site because of certain frustration. For discussions on 
square lattices with odd numbers of atoms per site, we 
will take this as a starting point and analyze disordered 
phases in section VIII B3. 

In the strong coupling limit where the hopping inte- 
gral is much less than Eg.c, after integrating over 0(2) 
and 0(3) degrees of freedom, or chargons and spinous, 
one obtains an effective low energy theory. 



□ □ □ □ 



(80) 



a^i is an Ising field defined at each link kl; □ is an el- 
ementary plaquette in a (1 -|- d)D Euclidean space. 
and are carried over spatial plaquettes and plaque- 
ttes involving temporal links respectively, fc,.,- can be 
estimated using a high temperature expansion 



ks (^^) 5 



(81) 



is a high energy cut-off in our problem. 
In one- dimensional lattices, the term proportional to 
ks is absent; furthermore, the result in Eq.80 is valid as 

far as 

t < GccEc- 

And therefore conclusions on one-dimensional lattices ar- 
rived in this section are expected to be valid for all one- 
dimensional Mott states disregarding the ratio between t 
and Eg. 

It is worth emphasising that only the topological term 
explicitly depends on numbers of particles at each site. 
General effects of Berry's phases on spin liquids have been 
studied in various papers on strongly correlated electrons 
[53,26-30,33]. The topological term in Eq.80 indeed has 
rather surprising consequencies on many-body ground 
state degeneracies and quantum numbers of excitations. 
In fact, this term determines the confining-deconfining 
property of effective Ising gauge fields in one-dimensional 
lattices. It can be easily appreciated in terms of suppres- 
sion of instantons. 

Consider a .^2-iiistanton in (1 -|- l)d Euclidean space 
which emits a string of negative bonds pierced by the x- 
axis so that the boundary condition along the temporal 
direction is periodic; only at the Xth plaquette shown in 
Fig.8, 



(82) 



otherwise the plaquette integral over an elementary 
square is one. 
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FIG. 8. A .Z2-instanton located at Xth plaquette in 
(1 -|- 1)D. At links pierced by the dashed line, cr?j = —1. 

Following the action in Eq.80, an instanton located at 
Xth plaquette along the a;-axis has a finite action and 
carries a phase 



Tb = XnNo. 



(83) 



In a dilute instanton gas limit (o( e '^^)), the partition 
function of instantons is 



n ■ {X^} {X„} 
n 

m=l, 



(84) 



where the subscript m labels instantons. For = 2n+l, 
the sign of the amplitude of an instanton therefore alter- 
nates between positive and negative ones as the location 
of the instanton is shifted by one lattice constant a along 
the X-axis. A direct evaluation shows that in this case in 
the leading order of o(e~'^^). 



InZ 



DIG 



PL, 







(85) 



as the perimeters along x and r direction in the Euclidean 
space go to infinity. It suggests that instantons of Z2 type 
should be absent in (1 -|- 1)D because of destructive in- 
terferences [61]. 

The alternating sign of amplitude also implies that the 
ground state should break the crystal translational sym- 
metry. Furthermore, the amplitude is periodical as a 
function of centers of instantons with a period of 2a. So 
the ground state should develop a spin-Peirls structure 
with a period doubling that of underlying lattices. The 
suppression of instantons therefore appears to lead to 
unusual deconfining Ising gauge fields in (1 -|- 1)D and 
a twofold degenerate spin-Peirls state, to which DVBCs 
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belong. For an even number Nq, obviously instantons 
proliferate in the system and the ground state is non- 
degenerate and translationally invariant. 

The instantons' behavior in this limit also determines 
quantum numbers of low lying elementary excitations, 
especially "topological" ones. For an odd number of 
bosons per site, a kink shown in Fig. 4 should be an ele- 
mentary spin-one excitation in the spectrum; this is in- 
dicated by the deconfining property of (frustrated) Ising 
gauge fields in (1+1)D. Thus, for an odd number of atoms 
per site, spin-charge separation takes place in any one- 
dimensional Mott state. On the other hand, for even 
numbers of atoms per site, gauge fields are confining. 
Lowest lying elementary excitations carry S = 1, Q = 1, 
or Q = 0,S = 2 or S = 0, Q = 2. The first one repre- 
sents adding one atom to the system; the second type of 
excitations is to flip one spin singlet into S = 2 state and 
the last one is to add (or remove) a singlet pair of two 
spin one bosons (Fig. 4c). 

Generally, topological terms are known to play very 
important roles in quantum disordered states. In one- 
dimensional Heisenberg antiferromagnetic spin chains, 
topological terms lead to gapless spin liquids in half- 
integer spin chains, differing from an integer spin chain 
[53]. In 2D, topological terms can result in spin-Peirls 
states with different periodicity for half integer, odd 
and even integer Heisenberg antiferromagnets (disor- 
dered states) [26,27]. In the current situation the topo- 
logical terms determine distinct quantum numbers of ex- 
citations for even or odd numbers of atoms per site. 

B. Quantum dimer models 

The Ising gauge fields describe the dynamics of low ly- 
ing collective spin states of condensates at each site. Let 
us consider 

= 4>\j't>io - \, 

4j = ^iK-Ai + KA-i - hl.ohio) (86) 

where (f)\j is the creation operator of a singlet state of 
two adjacent condensates, or a valence bond; hj^^ is a 
creation operator of a condensate with total spin one. 

In a projected space of valence bond configurations 
where each link is occupied by at most one valence bond 
defined by the Hilbert space at each link is an Ising 
doublet as pointed out in section V. One can easily con- 
firm that Tifj, a = x,y,z satisfy the algebra of trg and 
can be identified as erg-. So one interprets the effective 
action in terms of dynamics of valence bonds. 



An interesting alternative but closely related point of 
view is the quantum dimer model [32]. The effective 
Hamiltonian of the Ising gauge fields in Eq.80 can be 
written as 

{«> □ n 

where the first sum is over all links and the second sum is 
over spatial plaquettes; and in one-dimensional lattices, 
the second term is absent. One can verify that the action 
of HiG is precisely that of the Ising gauge fields in Eq.80 
[58]; indeed, 

ks = eKs,kr = In tanh(er) . (88) 

When r is much smaller than Kg, kr becomes much 
larger than unity. On the other hand, kr is much less 
than the unity as F becomes much larger than Kg. 

To generate desired topological terms in the action, we 
require that the Hilbert space satisfy the following local 
invariance at each site 

gh^ = ^,gk = exp[z^(^ + No)], (89) 

-I- 

for an even (A^o = 2n) or odd number (A^o = 2n + 1) of 
atoms each site. This constraint is identical to the con- 
straint in Eq.48 which was derived from a very different 
consideration. The sum -|- is carried over all links con- 
nected with site k. One can easily show that again for 
each site k 

gk'^niGgk = niG. (90) 

or TiiG is locally invariant under the action of gk- The 
parity of the number of bosons enters the theory only in 
the definition of the Hilbert space. 

At each link lives a two-dimensional Hilbert space 
cr^ = ±1. By assuming each link either occupied by a 
dimer (cr^ = ~1) or empty (cr^ = 1), one introduces 



as a dimer counting operator (see also section V). The 
first term in Hig can be interpreted as the chemical 
potential of dimers. As is a linear combination of 
creation and annihilation operators of dimers, in square 
lattices the last term is the "kinetic" energy of dimers, 
which conserves the parity of numbers of dimers emitted 
from each site but does not conserve numbers of dimers. 
The connection between Ising gauge fields and quantum 
dimer model was recently illustrated in an interesting 
work [33]. 

The Hilbert space has four distinct sectors if periodical 
boundaries are imposed along both x and y-directions. 
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This aspect of the gcnerahzcd model can be demon- 
strated exphcitly in a torus or ring geometry. Let us 
introduce winding number operators 

T.,y^l[(<yl,fly = l. (92) 

Here C^^^ is a path extending from — c» to +00 along 
the a:;(y)-dircction: the product is carried over all ver- 
tical links (y-direction) or horizontal links (x-direction) 
pierced by the path. Similar global operators were previ- 
ously introduced for the study of correlated states [31,33]. 
It is easy to verify that for i = x,y, 

[fi,fj]=0,[fi,niG]=0. (93) 

Tj = ±1 are good quantum numbers and the Hilbert 
space therefore has four sectors with different parities 
defined by T^^y. For a cylinder or a ring, one can only 
define Ty if x is the direction of circumference. In this 
case, there are two sectors in the Hilbert space. 

1 . Even numbers of atoms per site in all bipartite lattices 

Let us now turn to the ground staters of the Haniilto- 
nian in Eq.87. Solutions depend on the parity of numbers 
of bosons per site and we start with an even number of 
bosons per site. Following Eq.81, as t goes to zero, fcg.r 
approach zero; therefore in this limit, 

» 1. (94) 

The relevant Hamiltonian is 

W«d„^ = 2^^d; (95) 

which commutes with the dimcr counting operator de- 
fined at each link rj. The Hamiltonian is positive defined 
and the energy has a lower bound of zero. The ground 
state |G"^) obviously is a vacuum of dimers with d = 
at each link and is nondegenerate. In one-dimensional 
lattices, the ground state is 

d„\G^) = 0,fy\G^) = \G^). (96) 

where t] is an arbitrary link. 

In a torus, it has an even parity of T^^y; that is 

dr,\G^)=0X,y\G^) = \G')- (97) 

We identify this state as an SSMI state obtained in sec- 
tion II. The lowest excitations in this sector have energies 
4r higher than the ground state. One can also show that 
the lowest energy states in Ta; = — 1 or Ty = — 1 sector 
have energies proportional to the system size; they are 
not degenerate with the global ground state. 



2. Odd numbers of atoms per site in ID lattices 

When t is much less than Eg^c for an odd number A'o, 
the low energy Hilbert space becomes spanned by con- 
figurations where only one of all links connected with 
each site is occupied by a dimer. In ID, the second ki- 
netic term in Eq.87 is absent and the constraint leads 
to a dimcr crystal ground state with a twofold degener- 
acy. Particularly, two ground states [GJ 2) have different 
Tj^-parities 

dv\Gl2) = ^^^^^\Gl,),fy\Gl,) = ±\Gl^). (98) 

Here numerates link 77 in one-dimensional lattices. 
Interpreting each dimer as a valence bond, this point of 
view precisely leads to a dimerized-valence-bond crystal 
state. 



3. Odd numbers of atoms per site in square lattices 

For quantum spin disordered states in square lattices, 
one arrives at the following conclusions. For an odd num- 
ber of atoms, the number of dimers from each site should 
be 1,3,5... etc; however, as F goes to infinity, the low 
energy Hilbert space is spanned by configurations with 
one dimer emitted from each site. Consequently, the con- 
straint in this limit conserves the number of dimers emit- 
ted from each site and results in the following identity for 
any U{1) gauge choice 4>, 

hk = exp[i0(^ 4 - 1)] , /ifc* = (99) 

-I- 

The mapping between a number conserved dimer model 
and a compact C/(l) theory was obtained in [45]. 

In square lattices, the constraint for Nq = 2n+l atoms 
per site therefore still leaves an exponentially large de- 
generacy in the low energy manifold; correspondingly, the 
kinetic energy in the Hamiltonian becomes 

Kdm = -Ke (T^a-a^cj- + h.c. (100) 
□ 

in the reduced Hilbert space. This term lifts the residual 
degeneracy in the low energy manifold (cr^ are the rais- 
ing and lowering operators defined in cr^ basis.); on the 
other hand, it preserves the reduced Hilbert space 

K^'^ldm^k = '^Idm- (101) 

Terms such as crVW^, or a~a'^ are not invariant 

under the local U{\) gauge transformation and have zero 
matrix elements in the space defined in Eq. 99. 



16 



It can be shown that the Hamiltonian acting on a re- 
duced Hilbert space defined in Eq.99 is equivalent to the 
Rokhsar-Kivelson's quantum dimer model with V — 
[33,32]. The ground states in a square lattice at ^ = are 
column states with a fourfold degeneracy; furthermore, 
they break the crystal translational symmetry. This im- 
plies that a quantum disordered state in a square lattice 
should be a dimerized valence bond crystal state. For 
discussions on excitations, see appendix A. 

To summarize, two main approaches have been em- 
ployed to study spin correlated Mott insulating physics. 
The first one is using the effective Hamiltonian in re- 
duced Hilbert spaces. The problem of interacting spin- 
one bosons is first mapped into a constrained quantum 
rotor model; spin correlations in Mott states then are 
studied using the Hamiltonian Ti.Mi introduced in section 
II. In the strong coupling limit when J^x is much less than 
Eg, for an odd number of atoms per site, the Hamiltonian 
is further equivalent to the bilinear-biquadratic Hamilto- 
nian Hb.b. ■ 

The second approach is to map the quantum problem 
into a classical one; the CQR problem can be mapped 
into a classical 0{2) (g) 0(3) (8) Z2 model. In Mott states, 
the corresponding classical model is an 0(3) (g) Z2 one, 
with a topological term. At the strong coupling limit, the 
model is also equivalent to constrained quantum dimer 
models. This mapping is employed to investigate dis- 
ordered states in both one- dimensional and square lat- 
tices. For an odd number of particles per site, the later 
approach is used to generalize the results derived in the 
first approach to an intermediate hopping limit where the 
first approach is invalid. 

At last, in square lattices, the effect of topological 
terms when Nq is odd can be conveniently studied by 
a duality transformation from a constrained Ising gauge 
model to a fully frustrated Ising model in a transverse 
field [29,30,33,49]. The duality transformation indicates 
the ground state at infinit F limit should break the crys- 
tal translational symmetry. This again is consistent with 
the point of view of quantum dimer models. 

IX. EVIDENCE FOR SSQCS IN 
ONE-DIMENSIONAL OPTICAL LATTICES 

Following the analysises in section VI, it is clear that 
to have SSQCs in high dimensional lattices for integer 
numbers of particles per site. Eg has to be at least com- 
parable to Ec- On the other hand, following analysises 
in [34], because of long wave length fiuctuations, in one- 
dimensional lattices spin correlations are always short 
ranged; the rotational symmetry is unbroken (see Fig. 9). 
Spin excitations are fully gapped for any finite value of 
Eg. For these reasons, 

Gsc = +00. (102) 



For any hopping t and any Eg, Gg is always less than 
Gsc- Following discussions of Intermediate Hopping 
limit II in section VII, we therefore expect that when 
Gc > Gcc or t > EcGcc, ground states must be conden- 
sates of fractionalized chargons, either charge-e SSQC 
or charge-2e SSQC. These condensates are spin singlets, 
differing from conventional pBECs discussed in high di- 
mensional lattices (See Fig. 9 and Fig. 10). 

The effective Hamiltonian for one-dimensional lattices 
can be written as 

'^}U = -iT^^blbial^ + h.c.) 

(ki) 

+ T.^-^'P>' + ^^T.''M (103) 

k (kl) 

which differs from Ti-fqcb m Eq.66 by a term involving 
spatial plaquettes. The Hilbert space of this Hamilto- 
nian is subject to the same constraint as Eq.68 

Ci'^ = Ci' = exp {^.[blbu + ^ ^^]), 

+ 

and {blhk) = No. 

SSMI I SSQC(2e) 
(a) t/Ec 

DVBC I SSQC(e) 
'^'^ 

FIG. 9. States present in one-dimensional lattices with an 
even or odd number of spin-one bosons per site as the hop- 
ping t is varied, a) is for an even number of bosons per site;, 
b) is for an odd number of bosons per site. DVBC states 
in one-dimensional Mott states are suggested by the quan- 
tum dimer model in section VIII; in the limit t ^ ^/EsEc, 
these states are also suggested by the bilinear-biquadratic 
S = 1 spin chain model. We believe these states appear 
in cold atoms in an optical lattice when laser intensities are 
varied. Note that SSQCs have only quasi-long range order 
in one-dimensional limit and charge-e or -2e character might 
only be distinct close to the critical points. 

Therefore, the physics of spin-one bosons with anti- 
ferromagnetic interactions in one-dimensional lattices is 
effectively equivalent to the physics of spinless bosons in- 
teracting with each other via Ising gauge fields. This is 
the key idea behind the notion of quantum condensates 
of spin-one atoms; in this limit, spin-one atoms behave 
as if they lost their identities as spinful particles. 

Let us start with an even number of particles per 
site. In this case, the only fixed point of the Ising gauge 
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field theory (unfrustrated) in one-dimensional lattices is 
a strong coupling fixed point regardless the value of r^. 
The ground state for Ising fields should be unique with 
Ty — 1 and for any link kl, 



1. 



(104) 



And all excitations of Ising fields are fully gapped. The 
effective Hamiltonian can be reduced to TLpb in Eq.73, 
i.e., 

Ki = T.(^WMi + h.c.) + E ^ - ^Z^'^' 

k 



(kl) 



(105) 



Therefore when t2 < Ec, the ground state is a spin sin- 
glet Mott state, or an SSMI; furthermore, we believe that 
when t2 > Ec, chargons are paired and "condense" form- 
ing a charge-2e SSQC. This leads to a phase diagram 
for one-dimensional optical lattices in Fig. 9 a). 

For an odd number of particles per site and when 
Ec S> Es, one can carry out similar discussions. Fol- 
lowing discussions in section VIII A and B, we find that 
the ground state for Ising fields has a two-fold degeneracy 
with Ty = ±1, and for any link kl 



i) all = (-1)', or ii) all = ("1) 



fc+i 



(106) 



corresponding to two states breaking the crystal transla- 
tional symmetry. And excitations are gapped. In addi- 
tion, because of the Berry's phase term (see discussions 
in section III), the Ising fields are weakly interacting in 
this case. 

Hopping of a chargon in this gauge field background 
is similar to hopping of an atom in the projected spin 
singlet Hilbert space. A chargon hopps from site k to 
site fc -|- 2 via a virtual excitation involving an energy of 
Es- More specifically, in the representation introduced 
in section V 

one confirms that the hopping corresponds to the follow- 
ing process 

|...2n-H2,0,2n+l,l,2n-(-l,0...) 
|...2n + l,l,2n + 2, l,2n+l,0..> 
|...2n-|-l,l,2n-f l,0,2n-|-2,0..). 



At Gc > Gcc, chargons condense and the ground state 
is a charge-e SSQC, similar to condensates of charged 
solitons proposed for non-integer numbers of atoms per 
site [19]. Subject to strong one-dimensional fluctuations, 
SSQCs should be understood as states with quasi long- 
range phase order instead of condensates. 



(a) 





(b) 
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(c) 



FIG. 10. Condensation of particle-hole pairs suggested by 
solutions to Eq.l03. For an odd number of particles per site 
in a), a pair of particles with ±e charges are shown to con- 
dense; in b) for an even number of particles per site, a pair 
of particles with ±2e charges condense. Note that these con- 
densates are spin singlet states. As a reference, in c) we also 
show condensation of particle-hole pairs in a pBEC which is 
a stable phase only in high dimensional lattices. Locations of 
charges are indicated by dashed lines. 

The difference between condensates for odd and even 
numbers of atoms per site can be attributed to the selec- 
tion rule in low energy Hilbert spaces for each individual 
site defined in Eq.ll. In one-dimensional lattices, the 
selection rule defined in section II leads to a two-fold de- 
generate DVBC for an odd number of atoms per site and 
a non-degenerate SSMI for an even number of atoms per 
site in the strong coupling limit. Therefore, elementary 
charged excitations are either charge-e ones or charge-2e 
ones depending on the parity of numbers of atoms per 
site in one-dimensional lattices (see also Fig. 4). Conse- 
quently, condensates indicated in the fractionalized rep- 
resentation involve condensation of chargons or chargon- 
pairs depending on the parity of Nq. 



Consequently, the effective Hamiltonian should be, in- 
stead of a paired hopping form 



X. POSSIBLE BERRY'S PHASE EFFECTS ON 
SSQCS IN ONE DIMENSIONAL LATTICES 



Pk_ 

2C 



fipk- The action for Hamiltonian Tifq^b under the constraint 

defined after Eq.lOO can be derived following a standard 
(107) procedure; the result is 
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S = -Yl '^k,i<^ki cos(xfe - Xl) 

(kl) 

an k 

(108) 

Jk,i is given by 

Jk,k+x = {et), Jk,k+r = {eEc)-^ (109) 

and kr is given in Eq.81. This action should be used to 

study long wave length physics in one- dimensional lat- 
tices for an arbitrary hopping amplitude. 

When Jk,i is small, the phases Xk are disordered; these 
disordered phases represent Mott insulating states in one 
dimensional lattices which have been studied in details in 
section VIII. From a point of view of the action in Eq.108, 
differences between an even and odd number of atoms 
per site are only manifested in Berry's phase terms. The 
interesting "even-odd" effect in one-dimensional Mott 
states of spin-one bosons represents a well-known Berry's 
phase effect on spin liquids. 

Results in the previous section appear to imply that 
properties of "Higgs" phases close to critical points 
should also be solely determined by the even-odd parity 
of numbers of atoms per site. Especially, for A^o = 2n, 
the ordered phase represents a charge-2c SSQC and for 
No = 2n + 1, the ordered phase represents a charge- 
e SSQC. Whether an one-dimensional lattice with cold 
atoms perhaps is one of few systems where the Berry 
phase determines how the U{1) symmetry is broken spon- 
taneously remains to be understood. In a separate paper, 
wc arc going to numerically study SSQCs and the specu- 
lated effect of Berry's phases on spontaneous symmetry 
breaking. 



XI. CONCLUSIONS 

In this article, we have illustrated the notion of frac- 
tionalized atoms and demonstrated the possibility of hav- 
ing SSQCs of sodium atoms in one-dimensional optical 
lattices. Three major arguments have been developed 
to investigate spin singlet Mott states and condensates. 
These are: 

i) In a projected spin singlet Hilbert space, spin-one 
bosons with antiferromagnetic interactions are equivalent 
to spinless bosons interacting with Ising gauge fields. 

ii) Spin singlet Mott states are fully characterized by 
even- and odd-class quantum dimer models. 

iii) Because of a selection rule in the Hilbert space, in 
one-dimensional lattices, superfluid phases of atoms can 
have either charge-e or charge-2e characters depending 
on the parity of numbers of bosons per lattice site. 



It is worth emphasizing that charge-e and charge-2e 
condensates can be distinguished by studying frequen- 
cies in the ac Josephson effects. Finally, the even-odd 

parity effect discussed here also appears to be attributed 
to Berry's phase effects on spontaneous symmetry break- 
ing. 

Finally, we would like to mention two open questions 
which we believe deserve further investigation. First, 
the nature of states for non-integer numbers of parti- 
cles per site is not clear to us; possibilities of having 
spin singlet condensates in certain limits in this case ex- 
ist but need to be clarified. There are also issues con- 
nected with phase separation, similar to what happens 
in doped antiferromagnets [62,63]. Second, it is appeal- 
ing to verify results proposed in this article in the con- 
text of Luttinger- liquid theories; at the moment we are 
not aware of such attempts and do not know if results 
discussed in this paper can be rederived in a Luttinger- 
liquid theory based approach. Although it is not clear to 
us whether the phenomena discussed in this article are 
related to the Luttinger- liquid physics, attempts along 
this line of thought might shed light on the further un- 
derstanding of one-dimensional phases of spin-one bosons 
with antiferromagnetic interactions. 
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APPENDIX A: CHARGE-2E EXCITATIONS IN 
DVBCS IN SQUARE LATTICES 

For a DVBC in square lattices, the wave function can 
also be represented by Eq.43; the product is however car- 
ried over all even (odd) links in each row or in each col- 
umn. The ground state thus has a fourfold degeneracy 
[19]. 

The excitations in 2d DVBC states have some fascinat- 
ing properties. It is convenient to employ the following 
modified creation and annihilation operators of valence 
bonds. 
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where 0^ , 0^ are the creation and annihilation operators 
of spin singlet states of two coupled condensates defined 
in Eq.86 at a link 77 and 



^ , 1, if 4>14>.,^ - VP; 

'0, otherwise 



(Al) where is a creation operator defined at site k. M. is 
a path starting at site k and terminated at infinity and 
rjo^e are the odd and even links along the path (counted 
from site k). In one-dimensional lattices, this precisely 
creates a spin-one kink as shown in [19]. 

In a leading order of r/|y|, we estimate the energy of 
the spin-one excitation 



(A2) 



Finally, 5{r}) — 1 — 5{rf). The reduced Hamiltonian can 
be written as 



'vi (I'm (i'v2~T( 4>l^ 4)l^_ (f)fj^ (f)fj^ + h.c. ) ] 



(A3) 



771, 772 are a pair of two parallel adjacent links, perpendic- 
ular to another pair fji, 772 in an elementary plaquette. 

To facilitate discussions, we introduce a local gauge 
transformation Ek at site k = (m,7i); 



Ek = exp (i7r[^ - 1]) 



(A4) 



E : 



VLm 



(A8) 



where L^vi scales as the size of the system and the exci- 
tation is infinitely massive in a thermodynamical limit. 



c 



where the sum is again over all neighboring sites of site k. 
It is easy to confirm that all the resonating valence bond 
configurations and the Hamiltonian are invariant under 
this local gauge transformation 



rd- 



(A5) 



To facilitate discussions, we consider a limit where two 
parallel valence bonds attract each other strongly and 
—V is much greater than the exchange interaction T. 
(The following discussions are also expected to be valid 
as far as — < T because a bare — ^ is always renormal- 
ized to infinity in this limit.) By making an expansion 
over T/jT^I, one anticipates the following ground state 
wave function. 



\G1)^\CS^) + J2{^f\CS,) 



M 



V = -T ^ 4>l,^l^(i>fn(f>fi2 + h.c. 



(A6) 



where \CSi) is one of the four fold degenerate column 
states(see discussions at the beginning of this section). 

Let us start with solitonic spin-one excitations. Such 
a spin-one excitation is created by breaking one valence 
bond and sending one of the unpaired condensates to in- 
finity while keeping the other one fixed. The many-body 
wave function of a solitonic excitation localized at site k 
would thus be created by the operator 



n ' 



n ' 



(A7) 



FIG. 11. Moving an S = 1, Q = excitation (a dot with a 
double headed arrow) around a dual particle (defined along 
path C, see Eq. All) results in a minus sign in the many-body 
wave function; a dual particle and an S = 1, Q = excitation 
view each other as half vortices. 

The absence of spin-one solitonic excitations in the low 
energy excitation spectra doesn't depend on the partic- 
ular form of wave functions we used. It is a topologi- 
cal property of the ground state which determines the 
interactions between solitonic excitations. Here we are 
going to provide some understanding based on the du- 
ality between spin excitations which are "electric" and 
topological excitations with dual "magnetic" charges. 

To identify this gauge transformation Ek as the one of 
Ising gauge fields, one notices typically 



(A9) 



Here 7r(77, fc) is unity if link rj is connected with site fc; 
otherwise it is zero. As satisfies the algebras of Pauli 
matrices and should be identified as of the Ising gauge 
fields, the local gauge transformation effectively yields 
(T^ — CT^ and generates a desired Ising gauge transfor- 
mation. 

An operator carries an " electric" charge if it transforms 
nontrivially under the local gauge transformation. In our 
case, we can show that 



E^'^^^{k')h 



-1) 



(AlO) 



and the creation operator of a spin one solitonic exci- 
tation carries a charge defined with respect to the local 
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gauge transformation. The corresponding particle thus is 
topological in nature and violates the invariance defined 
in Eq.A5 only at the site k' where the particle is created. 

Furthermore, wc can introduce a creation operator of 
a dual magnetic charge; 

n^{x,y) = ]\{24>\4>r,-l) 

(All) 

where the path C begins at the center of a plaquette (a;, y) 
and ends at infinity by piercing all links intersecting with 
it. In the leading order of r/|y|, one also finds that 



vl/^+|G?>«.(-l)^-|G?) 



(A12) 



where Nm is the number of unoccupied links along the 
path A4. This implies that the ground state should be a 
condensate of dual particles. 

Now one can move a spin-one solitonic spin excitation 
^^j^ around a dual magnetic charge \E'^. In the Hilbert 
space spanned by valence bond configurations, the fol- 
lowing commutation relation holds: 

{{24,l4>r, - 1),4>1,} = {{24>l,4>r, - l),4>r,} = 0. (A13) 

Following Eq. A13, one verifies that the many-body 
wave function changes sign after this operation (shown 
in FIG. 11). This implies remarkable long range correla- 
tions in the system. Since the "magnetic" and "electric" 
charges in the theory see each other as half-flux-quanta 
, it is imaginable that as one of the charges condenses, 
the other one gets confined following a standard picture 
in topological field theories. In the current context, the 
dual magnetic charges condense and spin-one solitonic 
spin excitations are confined because of random phases 
developed when moving in a dual magnetic charge con- 
densate. 

We are now ready to address the quantum numbers of 
elementary spin excitations. Because of confining forces 
between S ~ 1 solitonic excitations, only bound states 
of soliton-anti soliton pairs can exist in the excitation 
spectrum as elementary excitations. So there will be two 
branches spin excitations with S = 1 and S = 2 respec- 
tively. Taking into account Eq.33, we should further ex- 
pect the bound states of S" = 2 have lower energies than 
S = 1. In fact, microscopically these bound states are 
likely to be the excited states of S' = 2 and S = 1 oi two 
coupled neighboring condensates. 

A charge-e soliton can be created by adding an atom 
to a site while keeping the rest of sites paired; this excita- 
tion is a spin singlet as the ground state but with an extra 
atom or charge-e. For the same reason discussed before 
for a spin-one solitonic excitation, a charged soliton is 
infinitely massive because of the long range rearrange- 
ment of valence bonds and intends to pair with another 
charged soliton. Naturally, at low energies, one can only 
have charge-2e spin singlet excitations in this case. 



APPENDIX B: THE DERIVATION OF THE 
FRACTIONALIZED REPRESENTATION 

The equivalence between Eq. 6 and Eq.56 can be 
demonstrated by examining the partition function, 



Z = Tr[exp{-(3nFR)V]; 



(Bl) 



■p is a projection operator defined in Eq.(B6) 

To facilitate the calculation, we introduce the following 
special "coherent" state representation 



No+SN 



(n • at)" 



No+SN 



{zb^y 



|o), 



0) 



(B2) 



assuming that A^o ^ but both are much larger than 
unity and z = exp{i(j)). One can easily show that at large 
N limit, 



{n\n')^6{n-n'),{z\z')^S{z-zy, 



(B3) 



and 



P\z) 



dlnz 



Sin) = inx -r— In). 

an 



(B4) 



For the study of the action, we slice the (d+l) dimen- 
sional Euclidean space into M slices along the temporal 
direction and rewrite the partition function as 

r 

\{nk,T+e}, {Zk,T+e}, Wk,r+t})- 

(B5) 

Here r and t + e label two adjacent slices and e = /3/M. 

The projection operator is to impose a constraint due to 
the symmetry of many-boson wave functions. 



P = ^ exp[i — ^7r(nfe + Sk)] 



(B6) 



so that only states with an even rifc -|- Sk contribute to 
the functional integral. 

The matrix element in Eq.B5 can be conveniently eval- 
uated by inserting a complete set of nk,Sk, the eigen 
states of the number and spin operators. 

^ {{nk,r}, {zkA^ Wk,T}\ exp{-e^FR)r ® 

{nfc},{-Sfc} 

\{nk}, {Sk}){{nk}, {Sk}\{nk,r+eh Uk,r+e}, {a^,.+J). (B7) 
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Here Zk = cxp{ixk)- 

Consider a simplified situation of a planar n = 
(cos 0, sin ^,0) (a generalization appears to be possible). 
As suggested in Eq.B4, S and p arc conjugate variables 
of n and 4>- So we can express the the eigen states \nk), 
\lk) in the following simple forms 



27r 

(nfe|5fe) = — ^exp(i(?ife5fe). 



(B8) 



The rest of the derivation is identical to that in [11]. 
Redefining as along a link located at site k between 
slice r and r + e and summing up all contributions at dif- 
ferent slices, we obtain the result in the Euclidean space 

as 

S = -^ Jrr'CTrr'COSXrr' " ^ Jrr' C0s(2Xrr') 

rr' rr' 



(B9) 



Here summation goes over sites r = (A;,t) in the space- 
time lattice, Xrr' — Xr — Xr', ^nd 

Qf = - 5"V3 

is a ncmatic order parameter, and a^,., = ±1 is an Ising 
field that lives on the links and rr = r + e. 
The coupling constants are 



1 



1 



JC TS tc,s c+M 

•^r,r±f — ^ g> "^r.ribf ~ ^ g' r,r±{x,v,z} ~ tt-iVQ, 

-etNo/4: 



j2c _ t2s _ n t2c,2s 



This is precisely the action of Eq. 1 derived in an early 

work [11] (the Berry's phase term was omitted there). 
Therefore we established an anticipated equivalence be- 
tween Eq. 1 and Eq.56. Finally, we would like to mention 
that Ising-gauge-thcory based approaches have also been 
employed to study electron fractionalization in previous 
works [28-31,33,47-49]. 



APPENDIX C: NMIS IN HIGH DIMENSIONAL 
LATTICES: REVISIT I 

As chargons are fully gapped, states in this limit are 
incompressible. Integrating out chargons' or 6-particles' 
degree of freedom in Eq.62 we end up with an effective 
Hamiltonian 



W/9ca = -iJ2i4a^la(rii + h.C.) + ^ ^ 



(kl) 



(kl) □ □ 



(CI) 



Calculations indicate that Fq, K^a should be functions of 
t, Ec- When t is much smaller than Ec, or Gc = tE~^ is 
small, the ratio KsaJ^a^ is much less than unity, or 



^ « i(ra 



Ec). 



(C2) 



However, we speculate that close to the critical value Gcc, 
the ratio in Eq. C2 should be divergent (one should also 
expect terms involving large loops) . Finally, the integra- 
tion also leads to a new constraint on the Hilbert space, 

Cfa^ = Ct = exp (m[at a, + J2 • (C3) 



The sum over " -|-" is again carried over all links connected 
to site k. And the Hamiltonian is locally invariant under 
the action of (7^" 



fa 



n 



fqca- 



(C4) 



The condensation can be easily visualized when t is 
close to EcGcc, or in the large Ksa limit. The gauge 
fields are weakly interacting and fluctuations are small. 
For instance, one can estimate the probability of finding a 
loop of plaquettes with Yl^ a^.i = — 1 in (2-|-l) Euclidean 
space is exponentially small. 

It is therefore tempting to work out a self-consistent so- 
lution in a "mean field" approximation where Ho ""w = ^ 
at each plaquette. Furthermore, one chooses a unity 
gauge 



all = 1 7 fo'' liiik. 



(C5) 



One then obtains a FQCa, or an a-type FQC, as the 
ground state in this case. In FQCa, 



g''{k, fc + oo) = 0, Gafiik, k + oo) n„n/3 



(C6) 



as spin-one spinous condense in a state created by an 
operator aj^^g • n. 

It is also possible to show the existence of the fraction- 
alized phase beyond this simple mean field theory. In a 
dilute gas approximation suggested in [57], one finds an 
up-bound of the correction to G^isi^: k + oo) calculated in 
a mean field approximation which is exponentially small 
as Ksa or Gg is large. One therefore expects that the 
conclusions arrived above are valid beyond the mean field 
approximation. 

The ground state of FQCa does not exhibit phase co- 
herence or Josephson effects. However, the rotational 
symmetry is broken. There are two branches gapless spin 
wave excitations (spin one); the zeroth sound involving 



22 



the compression or expansion of particle-densities is fully 
gapped. 

A more tricky situation is when t is much less than Ec 
or correspondingly when Ta is large and gauge fluctua- 
tions are strong in the absence of matter fields. To argue 
the existence of condensation in this case, one recalls that 
there arc no phase transitions when Ka is varied but Gs is 
kept at infinity. Physically, presence of matters strongly 
renormalizes gauge field dynamics. A condensate repre- 
sents a self-consistent solution to the problem provided t 
is much larger than Es ; this conventional point of view is 
supported by results in [64]. However, since pure gauge 
fields at large Ta limit are confining, only pairs of a^a^ 
are free excitations at small t limit; as t is increased, con- 
densation appears to only involve these pairs [65,66]. To 
illustrate this point, we present some detailed discussions 
on pFQCa. 

As stated at the beginning of this section, there are 
even numbers of atoms at each site. Consequently, each 
site is occupied by an even number of chargons; excita- 
tions involving a change of numbers of chargons at each 
site are gapped by an energy of order of Ec. At ener- 
gies much lower than E^ the constraint in Eq.55 then 
indicates that the low energy Hilbert space should be 
spanned by states satisfying 



*,C^ = exp(i7ra[-afe). 



(C7) 



In this subspacc defined by two-particle operators 
(aj^.)^, the relevant Hamiltonian is due to pair hopping 
following Eqs.56,62, 
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(kl) 

(Cfc) T^paCk = '^pa- (C8) 



In a pFQCa, I) vanish as k ~ I ap- 

proaches infinity. However, the rotational symmetry is 
still broken and there exists the following long range or- 
der in the generalized BCS wave function 



A«^(fc,fc) = 3Q„, 



(Cll) 



which is invariant under the local gauge transformation 
defined in Eq.58. 

The diff'erence between pFQCa and FQCa is subtle. 
Both of them break rotational symmetry, support gap- 
less modes and have similar local dynamics. On the other 
hand, topological excitations in two states are distinct. 

Let us also mention briefly that when Eg ^ t,ti, only 
localized singlet pairs of spinous are allowed in the sys- 
tem. These states can be identified as SSMIs discussed 
in section II. 



APPENDIX D: DVBC'S IN A FRACTIONALIZED 
REPRESENTATION 

As shown in some details, the ground state for an odd 
A^o when Gs,t <C 1 is a projected superposition of paired 
condensates in {q, —q) and {q, it — q) channels and has a 
two-fold degeneracy. The corresponding wave functions 
of the two-fold degenerate ground states are topologically 
identical to 



^G.G.exp[^ 57C(77r)]|0); 



7=0,1 



and 



ti 



J^iTa^Eo). 

J- n. 



(C9) 



By examining the Hamiltonian in this reduced space, as- 
suming ti = iiN^ » Es, one concludes that the ground 
state is a condensate of paired spinous, or pFQCa. In a 
lattice with Vr sites, the wave function of a pFQCa can 
be written as (up to a normalization factor) 



i^i) = n 




^iv„xv.exp($t(Qo = 0))|0). 



(CIO) 



which explicitly satisfies the above constraint. Here 
^aiQo) is a creation operator for a pair of spinons in 
{Qo ~ Q, q) channel. The sum over q is carried over the 
first Brillouin zone. Vn^xVt is to project out Na x Vr- 
particle states. 



where h(q) is chosen to be cxp{iq). Two fold degener- 
ate states correspond to 50 = 5i = 1 or go = —gi = 1. 
It is straightforward to establish an equivalence between 
this projected pFQCa state and dimerized valence bond 
crystals studied previously [19]. 



APPENDIX E: THE DERIVATION OF EQ.80 

Let us consider a limit where x and n have short range 
correlations and 



{cOSXr COSXl) = C^, {-Orni) = Cn- 



(El) 



only if fc, / are two neighboring sites. 

Using the standard high temperature expansion tech- 
nique, one integrates out n,x in the action given in 
Eq.B9. The results in the leading order are 
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/ j.\4 \ z z z z 

-{et) 2^ <^kl(^lm<^mn<^nk 

kl,lm,mn,nk 

t t 



kl,l 



ra.mn.nK 



exp(ziVo^^^7r). 



(E2) 

(E3) 
(E4) 



The sum in Eq.E2 is over all elementary spatial plaque- 
ttcs occupied by links kl,lm,mn,nk: the sum in Eq.E3 
is over elementary plaquettes involved both spatial links 
kl, mn and temporal links Im, nk; without losing gc;ucr- 
ality, wc have sot C^^ and Cn to be unity. The structure 
of Eq.E2 is identical to Eq.80; further more, comparing 
Eqs.E2,E3,E4 with Eq.80, one obtains results in Eq.81. 



APPENDIX F: GAUGE INVARIANT 
CORRELATION FUNCTIONS 

It is tempting to classify SSQCs or more general FQCs 
in terms of local order parameters similar to those in 
Eq.l. To characterize FQCs, in addition one should in- 
troduce 



0° = i^kai^ka). 



In terms of order parameters C°'^'^, 

pBEC: Oi^O,0^^^o,OVO; 
SSQC(2e): 01=0,0^0,0^^=0; 
pFQCa: = 0, Ol^ ^0,0° = 0. 



(Fl) 



(F2) 



For FQCa and SSQC(e) states, the order parameters are 
nonlocal in terms of tp^ and will be introduced in a subse- 
quent paper on fractionalized states for non-integer num- 
bers of bosons per site. 

We would like to make two more remarks on FQCs. 
First, Elitzur theorem claims that in the presence of lo- 
cal gauge invariance, operators such as a^, a.\, 5^ which 
transform nontrivially under the local gauge transforma- 
tion can not develop nonzero expectation values [67] . The 
existence of FQCs which so far is explicitly demonstrated 
in a fixed gauge at first sight appears to be at odds with 
the Elitzur theorem. The apparent paradox can be for- 
mally resolved if one carries out a similar calculation in 
terms of gauge invariant correlators and confirms long 
range order. 

Indeed, for discussions on quantum condensates, one 
can also study the following gauge invariant two-point 
correlation functions 



/4 



Chi 



G\Kl) = {hl'[{alh){X{al)-' 

Cki ct 



/4 



(F3) 



where the products are carried over links along path Cki 
and C^i- Cki is a path connecting k and I site and 
is a closed path of four times as long. These correlation 
functions characterize the condensates when the gauge 
fields are weakly interacting [68]. 

Gauge invariant correlators suggest that it should not 
be bare a- or 6-particles but some nonlocal "particles" 
which condense. One such candidate is an a or 6-particle 
dressed in a string of gauge fields 

^LiCk) = aL n 4iCk) = fel n < 



Ck 



Ck 



where Ck is a path starting at site k and terminated at 
infinity. One can easily confirm the following local gauge 
invariance of these operators 



icir'Bict=Bi 



t 

ka ' 



(F5) 



These particles can condense without violating the 
Elitzur theorem. On the other hand they carry exactly 
the same charge and spin as bare particles of a^ and 6''' . 
We speculate that condensation of a or 6-particles in a 
fixed gauge might also imply condensation of particles 
created by and . 

The second remark concerns the relation between 
charge-c SSQCs and previously discussed spin singlet 
paired condensates or charge-2e SSQCs [11]. Both con- 
densates are rotationally invariant spin singlets and both 
are phase coherent. However, from the point of view 
of correlations it is obvious that charge-e SSQC doesn't 
involve pairing of particles and differs from paired con- 
densates. An explicit construction of wave functions of 
charge-e SSQC is difficult but available. At the time of 
writing, we believe that FQCs represent condensation of 
some topological solitons which in short are called as a-, 
b- particles. Particularly, charge-e SSQCs should exist 
in lattices with both non-integer and integer numbers of 
atoms per site. 
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